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1.  Introduction 


Many  conventional  projectile  configurations  contain  internal  parts  that  slightly  move  in  flight  in 
some  way,  shape,  or  form.  Fuze  mechanisms  used  on  some  indirect  fire  ammunition  employ  a 
rotor  that  is  permitted  to  move  slightly  with  respect  to  the  main  projectile  body.  Submunitions 
deployed  from  indirect  fire  projectiles  are  keyed  into  place  inside  the  round;  however,  small 
relative  motion  between  parts  occurs.  These  configurations  can  experience  dynamic  instability 
typified  by  large  loss  in  range  and  large  spin  decay  (i).  Soper  (2)  evaluated  the  stability  of  a 
spinning  projectile  that  contains  a  cylindrical  mass  fitted  loosely  into  a  cylindrical  cavity.  The 
cylinder  is  constrained  to  spin  with  the  main  body  projectile.  It  is  shown  that  an  unstable  coning 
motion  exists  in  which  spin  decay  and  cone  angle  grow  proportional  to  friction  coefficient 
between  the  mass  and  cavity  and  the  maximum  cant  angle  between  the  mass  and  the  projectile. 
Using  a  similar  geometric  configuration,  Murphy  (3)  developed  a  quasi-linear  solution  for  a 
projectile  with  an  internal  moving  part.  It  provided  an  explanation  of  the  unusual  flight  behavior 
exhibited  by  four  projectiles  that  all  contained  parts  with  slight  relative  motion  between 
components.  Later,  D’Amico  (4)  performed  a  detailed  series  of  experiments  where  a  projectile 
with  a  loose  internal  part  was  driven  by  the  rotor  of  a  freely  gimbaled  gyroscope.  The  gyroscope 
yaw  history  and  the  orbital  motion  of  the  loose  part  were  measured  and  used  to  predict  the 
moment  and  resulting  yaw  growth  caused  by  the  loose  part.  Hodapp  (5)  expanded  the  work  of 
Soper  and  Murphy  (2,  3)  by  considering  a  projectile  configuration  with  a  partially  restrained 
internal  member  with  a  mass  center  offset.  Results  of  this  study  indicate  that  small  mass  center 
offset  of  the  partially  restrained  internal  member  can  reduce  the  instability  caused  by  the  loose 
internal  part. 

New  projectiles  configurations  have  emerged  that  contain  multiple  moving  parts  which  are 
fundamental  to  the  basic  design  and  operation  of  the  projectile.  The  gimbal  nose  projectile 
configuration  is  an  example  of  one  such  configuration.  It  consists  of  a  standard  projectile  shape 
with  a  nose  section  that  is  free  to  rotate  with  respect  to  the  main  body.  Goddard  (6)  originally 
conceived  this  device  for  aircraft  control,  and  later  Barrett  and  Stutts  (7)  considered  this 
mechanism  for  active  control  of  munitions.  Schmidt  and  Donovan  (8)  as  well  as  Costello  and 
Agarwalla  (9)  showed  that  dispersion  of  a  fin- stabilized  direct  fire  projectile  could  be  reduced  by 
more  than  50%  using  a  passive  gimbal  nose  to  significantly  reduce  aerodynamic  jump.  Another 
multiple  component  configuration  is  the  dual-spin  projectile,  which  consists  of  forward  and  aft 
sections  connected  through  a  bearing,  allowing  different  spin  rates  for  each  section.  Smith  et  al. 
(JO)  used  a  dual-spin  projectile  for  active  control  of  an  artillery  shell  by  mounting  canards  on  the 
forward  section  of  the  projectile.  The  forward  section  was  roll  stabilized  to  aid  the  functionality 
of  the  canards  while  the  aft  section  provided  spin  stability.  Costello  and  Peterson  (JJ)  developed 
a  linear  theory  for  dual-spin  projectiles  that  predicts  stability  for  this  configuration  while 
Burchett  et  al.  (J2)  predicted  swerve  of  a  dual-spin  projectile  caused  by  lateral  pulse  jets. 
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An  internal  rotating  disk  is  an  important  dynamic  component  of  some  new  projectile 
configurations.  For  example,  a  new  concept  for  generating  real-time  bomb  damage  information 
relies  on  releasing  a  relatively  small  sensor  projectile  that  is  tethered  to  a  parent  bomb.  As  the 
two  projectiles  separate,  a  reel  on  the  parent  munition  spins.  In  another  application,  active 
trajectory  control  is  achieved  by  controlling  the  spin  rate  of  the  external  projectile  body,  thus 
predictably  changing  the  aerodynamic  loads.  To  implement  this  control  concept  requires  an 
internal  rotating  disk.  In  these  cases,  weapon  system  designers  require  guidance  on  the  effect  of 
the  rotating  internal  part  as  well  as  guidance  on  how  to  optimally  configure  such  a  system.  The 
work  reported  here  sheds  light  on  these  matters  by  first  developing  a  projectile  linear  theory 
specific  to  projectiles  with  an  internal  rotating  disk.  This  is  followed  by  a  study  of  the  effect  of 
the  physical  parameters  such  as  the  placement,  orientation,  mass,  and  speed  of  the  rotating  disk 
on  the  epicyclic  modes  of  vibration. 


2.  Rotating  Internal  Part  Projectile  Dynamic  Model 


A  projectile  containing  an  axisymmetric  rotating  internal  part  that  spins  at  a  constant  rate,  is 
considered  as  shown  in  Figure  1 .  The  mathematical  model  describing  the  motion  of  the 
projectile  allows  for  three  translational  and  three  rotational  rigid-body  degrees  of  freedom.  In 
order  to  develop  the  dynamic  equations  of  motion  for  these  six  degrees  of  freedom,  three 
separate  reference  frames  are  used  as  shown  in  Figure  1 .  The  ground  surface  is  used  as  an 
inertial  reference  frame  with  Kj  positive  down.  A  body  frame  is  fixed  on  the  projectile  at  the 
mass  center  of  the  two-body  system  with  /^positive  out  the  nose  of  the  projectile.  Due  to  the 
axisymmetric  nature  of  the  rotating  internal  part,  the  mass  center  of  the  internal  part  is  fixed  with 
respect  to  the  body  frame.  The  disk  is  considered  to  have  a  known  constant  spin  rate,  and  its  axis 
of  spin  is  specified  in  the  body  frame  by  a  set  of  direction  cosine  elements. 

The  three  translational  degrees  of  freedom  are  the  three  components  of  the  two-body,  mass 
center  position  vector. 

+  (1) 

A  sequence  of  rotations  from  the  inertial  frame  to  the  projectile  frame  is  defined  by  a  set  of 
body-fixed  rotations  that  are  ordered  in  the  conventional  manner  as  shown  in  Figure  2.  The  three 
rotational  degrees  of  freedom  are  the  Euler  roll  angle  ( ),  pitch  angle  {0),  and  yaw  angle  {¥). 

In  the  normal  process  of  simplifying  the  equations  of  motion  using  projectile  linear  theory,  an 
intermediate  frame  is  utilized,  namely,  the  no-roll  frame,  defined  as  an  intermediate  frame  before 
roll  angle  rotation.  Figure  3  shows  the  relative  locations  of  the  projectile  and  disk  centers  of 
gravity  and  the  projectile  body  centers  of  pressure. 
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Figure  1.  Position  coordinates  schematic  of  a  rotating  internal  part  projectile. 


Figure  2.  Attitude  coordinates  schematic  of  a  rotating  internal  part  projectile. 
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Figure  3.  Rotating  internal  part  projectile  geometry. 


The  kinematic  differential  equations  define  6  of  the  12  dynamic  equations  needed  to  describe  the 
motion  of  the  states:  x,  y,  z,  (p,  6,  and  y/.  The  no-roll  variables,  u,  v,  w,  p,  q,  and  r  ,  are  chosen 

for  the  remaining  12  state  variables.  The  transformation  from  the  no-roll  frame  (N)  to  the  inertial 
frame  (7 )  is 


4 

1 

>  = 

* 

II 

^  N 

\ 

1 

0 

1. 

(2) 


while  the  transformation  from  the  projectile  body  frame,  (5),  to  the  no-roll  frame,  (AO,  is 


1 

o 

1 - 

o 

^  N 

K, 

>  =. 
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_ _ _ _ _ 1 

_ 

B 

II 

Jb 

(3) 


leading  to  the  transformation  from  the  body  frame  (B)  to  the  inertial  frame  (I)  described  by 


1 

^e^y/ 

S4,SeCy,-c^s^  c^SgC^+s^s^ 

'Ib 

'h' 

h 

>  = 

* 

Jb 

Jb 

1 

_ 1 

^(jPe 

Ks 

Kg 

(4) 


In  the  above  equations  and  the  equations  shown  below,  the  standard  shorthand  notation  for 
trigonometric  functions  is  used:  sin(«)  =  ,  cos(«)  =  ,  tan(a)  =  . 
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The  mass  center  velocity  vector  of  the  two-body  system  is  defined  in  each  of  the  reference 
frames  previously  discussed, 

Vp//  =  +  vJ^  +  wKj^  =  uig  +  vJg  +  wKg  =  xlj  +  yJj  +  zKj ,  (5) 

as  is  the  angular  velocity  vector  of  the  projectile  body, 

^Bl I  ~  N  N  ^  N  ~  B  B  B  ’ 

and  the  angular  velocity  vector  of  the  no-roll  reference  frame, 

^N/J  ~  N  N  ^^NB  • 

The  angular  velocity  vector  of  the  internal  rotating  part  with  respect  to  the  inertial  frame  is  found 
by  summing  the  angular  velocity  vector  of  the  projectile  body  with  respect  to  the  inertial  frame 
with  the  angular  velocity  vector  of  the  internal  rotating  part  with  respect  to  the  projectile  body 
frame: 

^D/I  ~  ^B/I  ^D/B  ’  (^) 

where  the  angular  velocity  vector  of  the  internal  rotating  part  with  respect  to  the  projectile  body 
frame  is 

g^d/p  ~  B  B  n^ClKg  .  (9) 

Applying  the  transformation  given  in  equation  2  to  the  mass  center  velocity  components 
expressed  in  the  no-roll  reference  frame  yields 


Equating  the  projectile  angular  velocity  vectors  described  using  no-roll  frame  components  and 
using  Euler  angle  rates  generates 

^  10  tg  {p 

<^[=0  1  0  (11) 
\j/  0  0  Ijcg  r 

The  kinetic  differential  equations  are  derived  by  separating  the  two-body  system  at  the  disk  axle 
connection  points  and  considering  the  reaction  forces  and  moments  associated  with  each 
individual  part.  A  constraint  force  Fg  and  a  constraint  moment  M ^ ,  applied  at  the  disk  center 
of  gravity,  couple  the  disk  and  projectile  bodies.  The  translational  kinetic  differential  equations 
for  each  body  are  given  by  equations  12  and  13: 
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Summing  equations  12  and  13  yields  the  expression  for  the  translational  dynamic  equation  of 
motion  for  the  two-body  system, 


where 


and 


^C^C/I  ~  ^  ■■ 


W,=Wp+W^- 


'dV, 


cn 


dt 


■*"  ^NII  ^  ^C/1 


(14) 


(15) 


(16) 

(17) 


The  constraint  force  is  obtained  by  subtracting  equation  13  from  equation  12: 


F,= 


r?  \ 


^Djl  ^/>// 


V  J\ 

The  acceleration  of  the  mass  center  of  the  disk,  dp,/j ,  can  be  expressed  in  terms  of  the 

acceleration  of  the  projectile  mass  center  by  using  the  formula  for  two  points  fixed  on  a  rigid 
body. 

^D/i  ~  ^p/i  ^p/i  ^  -F  cOp/j  X  (^cOpij  X  fp^pi) .  (19) 


nip^mp 
ynipi  +  nip  j 


(18) 


After  substituting  this  expression  into  equation  18,  the  constraint  force  is  expressed  in  the 
following  manner: 


F,= 
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-F  (Opu  X  {(Opu  X  Pp^p, )  + 
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py 
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p  J 


(20) 


The  rotational  kinetic  equations  of  motion  for  the  projectile  and  disk  bodies  are  given  by 


'  dH 


DU 


dt 


=  M 


R  ’ 


’dH 


p/I 


dt 


—  M  p  Pp^D  ^  Fp+  M  ^ . 


(21) 

(22) 


Summing  equations  21  and  22  eliminates  the  reaction  moment  and  forms  the  rotational  kinetic 
equation  for  the  two-body  system  expressed  in  the  body  reference  frame. 


(23) 


^  'P^D  ^  ^  R  ^  A  • 


dt 


dt 


Unlike  the  kinetic  translational  equation  of  motion,  it  is  easier  to  form  the  rotational  equation  of 
motion  in  the  body  frame  and  later  convert  it  to  the  no-roll  reference  frame  when  expressing  in 
component  form.  Therefore,  the  angular  momentum  derivatives  are  expressed  as 
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DU 
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dt 

'dH 


PH 


dt 

Uh 


+  COg/i  X  H , 


PH 


+  COg/j  X  H  pii  . 


dt  dt 

The  translational  dynamic  equation  given  in  equation  26  is  expressed  in  the  no-roll  frame. 


(24) 
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where  the  weight  force  and  aerodynamic  loads  expressed  in  the  no-roll  frame  are  written  as 


(27) 

(28) 


so  that 
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Zw 

where 


>  =  mcg< 

0 

Zw 

(29) 


(30) 


The  aerodynamic  forces  are  described  in  section  3. 

The  body  frame  components  of  the  rotational  dynamics  equation  of  motion  given  by  equation  23 
can  be  written  as 
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where 
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The  unit  vectors,  Hx,  riy,  and  n^,  are  the  set  of  direction  cosine  elements  that  specify  the  axis  of 
spin  of  the  disk  in  the  body  frame. 

To  be  consistent  with  projectile  linear  theory,  equation  31  is  converted  to  the  no-roll  frame 
through  multiplication  of  the  equation  by  Tj .  Also,  a  change  of  variables  is  introduced  from 

body  frame  angular  velocity  components  to  no-roll  frame  components.  With  these  conversions, 
equation  31  is  expressed  as 


q  1  ^RD )’ 


where 


^RZ)  -  Bj{lp  +T^Ip,Tp,  -m^Sppp,Sppp,)r^  , 
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,  (38) 


(39) 


(40) 


In  the  cross  product  operator  matrix,  ,  in  equation  38,  the  body  frame  angular  velocity 
components  p,  q,  and  r  are  replaced  by  p  ,  c^q  +  s^r  ,  and  -  s^q  +  c^r  ,  respectively. 

Equations  10,  11,  26,  and  36  provide  12  nonlinear  differential  equations  that  govern  atmospheric 
flight  of  a  projectile  equipped  with  an  axisymmetric  rotating  internal  component.  With  a  given 
set  of  initial  conditions,  these  equations  can  be  numerically  integrated  forward  in  time. 


3.  Aerodynamic  Forces  and  Moments 


The  equations  of  motion  previously  discussed  are  largely  driven  by  the  aerodynamic  forces  and 
moments  exerted  on  the  projectile  body.  The  aerodynamic  loads  consist  of  steady  aerodynamic 
forces  and  linear  Magnus  forces  and  are  formulated  separately. 


'x. 

'x^' 

Y, 

>  =  < 

>  +  < 

Im 

z. 

Zm 

(41) 


The  steady  aerodynamic  forces  act  at  the  center  of  pressure  of  the  projectile  body  and  are 
provided  in  equation  42: 
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(42) 


6^X0  ■*"  6^X2^ 

CyQ  +  Cyg^fi 
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^ZQ  ^ZBl^ 

The  Magnus  force  acts  at  the  Magnus  force  center  of  pressure,  which  is  different  from  the  center 
of  pressure  of  the  steady  aerodynamic  forces.  Figure  3  shows  the  relative  locations  of  the 
projectile  body  centers  of  pressure. 


0 

pDCyp^a 

2y 

pDCypy^fi 


2V 


(43) 


The  longitudinal  and  lateral  aerodynamic  angles  of  attack  used  in  equations  42  and  43  are  given 
in  equation  44: 

«  =  tan  '(vv/m)  =vv/m  =  tan  '(v /m)  =  v /m  ;  (44) 

q^  =  '^pD\U^+v^+w^).  (45) 


Aerodynamic  coefficients  in  equations  42  and  43  depend  on  the  local  Mach  number  at  the 
projectile  mass  center. 


The  externally  applied  moments  on  the  projectile  body  found  on  the  right-hand  side  of  the 
rotational  kinetic  equations  contain  contributions  from  steady  and  unsteady  aerodynamics  and 
Magnus  moments. 
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(46) 


The  steady  aerodynamic  moments  are  computed  for  the  projectile  body  with  a  cross  product 
between  the  steady  body  aerodynamic  force  vector  and  the  distance  vector  from  the  center  of 
gravity  to  the  center  of  pressure.  Magnus  moments  on  the  body  are  computed  in  a  similar  way, 
with  a  cross  product  between  the  Magnus  force  vector  and  the  distance  vector  from  the  center  of 
gravity  to  the  Magnus  center  of  pressure.  The  unsteady  body  aerodynamic  moments  provide  a 
damping  source  for  projectile  angular  motion  and  are  given  by  equation  47: 
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2V 

Air  density  is  computed  using  the  standard  atmosphere  (13). 

4.  Rotating  Internal  Part  Projectile  Linear  Theory 

The  preceding  equations  of  motion  in  their  current  state  are  highly  nonlinear  and  while  a 
solution,  given  an  initial  set  of  conditions,  may  be  obtained  numerically,  it  is  desirable  to  solve 
them  with  a  closed  form  solution  for  increased  understanding  of  the  dynamic  behavior.  Linear 
theory  for  symmetric  rigid  projectiles  introduces  a  series  of  assumptions  that  yield  a  refined  set 
of  linear  differential  equations  that  can  be  solved  in  closed  form.  These  equations  form  the  basis 
of  classic  projectile  stability  theory.  This  same  set  of  assumptions  can  be  used  to  establish  a 
linear  theory  for  projectiles  containing  an  axisymmetric  rotating  internal  part  in  atmospheric 
flight.  The  necessary  assumptions  are  as  follows: 


•  The  variable  is  changed  from  no-roll,  station  line  velocity  u,  to  total  velocity  V.  Equations 
48  and  49  relate  V  and  u  and  their  derivatives: 


V  =  -i-v^  -1- ; 

(48) 

V  =  ipu  -1-  vv  -1-  ww)/V  . 

(49) 

•  The  variables  are  changed  from  time,  t,  to  dimensionless  arc  length  s.  The  dimensionless 
arc  length,  as  defined  by  Murphy  (14)  is  given  in  equation  50  and  is  measured  in  units  of 
distance  traveled: 

1  ^ 

s  =  —  f  Vdr . 

D  ^ 

(50) 

Time  and  arc  length  derivatives  of  a  dummy  variable  ^  are  related  by 

C  =  (VID)C; 

(51) 

^  =  (v  iDf((;"  +  <;v'  IV). 

(52) 

•  Euler  yaw  and  pitch  angles  are  small. 

«  1  s^^y/  «  1 . 

(53) 
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Aerodynamic  angles  of  attack  are  small. 


a^wlV  p~vlV. 


(54) 


•  The  effect  of  mass  and  inertia  changes  of  a  projectile  on  stability  are  well  known.  To 
properly  evaluate  the  effects  on  stability  solely  due  to  the  rotating  internal  part,  the  total 
mass  of  the  two-component  system  is  held  constant  and  individual  component  mass  and 
intertia  properties  are  appropriately  modified.  This  is  accomplished  by  reevaluating  the 
inertial  properties  of  the  projectile  after  the  addition  of  the  internal  part  based  on  its 
nominal  characteristics  prior  to  the  addition  of  the  internal  part.  Equation  55  provides  the 
shift  in  the  projectile’s  center  of  gravity  due  to  the  addition  of  the  internal  part,  and 
equation  56  is  used  to  populate  the  projectile’s  inertia  matrix  for  various  disk 
configurations. 
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where 


nip  =  nip  -  nip,. 

(57) 

•  The  projectile  is  aerodynamic  ally  symmetric. 

C  =  C 

'^NR  ^MQ  ■ 

(58) 

C  =C  =0 

(59) 

c  =  c  =  c 

'^YB\  '^ZB\  '^NA  ■ 

(60) 

•  A  flat  fire  trajectory  is  assumed,  and  the  force  of  gravity  is 

neglected  for  stability  analysis. 

•  The  quantities  V  and  (j)  are  large  compared  to  6,  y/,  q,  r ,  v, 
small  quantities  and  their  derivatives  are  negligible. 

and  w ,  such  that  products  of 

After  application  of  these  assumptions,  the  differential  equations  can  be  solved  in  closed  form  as 
stated  previously. 

x'  =  D. 

(61) 

y'  =  {dIv)v  +\i/D  . 

(62) 
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Equations  61-69  are  linear,  except  for  the  total  velocity,  V,  which  is  retained  in  several  of  the 
equations.  It  is  assumed  that  V  changes  slowly  with  respect  to  the  other  state  variables,  and  is 
considered  to  be  constant  where  it  appears  in  other  dynamic  equations.  With  this  assumption, 
the  total  velocity,  the  angle  of  attack  dynamics,  and  the  roll  dynamics  all  become  uncoupled, 
linear-time  invariant  equations  of  motion. 

Equation  69  is  the  matrix  form  of  the  epicyclic  dynamic  equations.  The  equations  largely 
determine  the  angle-of-attack  stability  of  the  projectile.  It  should  be  noted  that  the  basic 
structure  of  the  epicyclic  dynamic  equations  of  a  rigid  projectile  are  the  same  as  previously 
shown.  Differences  in  the  epicyclic  dynamics  of  both  configurations  are  contained  in  the 
coefficients  B,  C,  E,  and  E.  In  the  most  general  case,  where  the  disk  is  located  off  the  axis  of 
symmetry  and  canted  at  an  arbitrary  angle,  these  coefficients  are  algebraically  lengthy.  While 
straightforward  and  computationally  trivial  to  compute,  space  limitations  here  prevent  the  most 
general  form  of  these  coefficients  to  be  listed.  However,  Appendix  A  provides  these 
coefficients,  as  they  would  appear  in  the  C -programming  language.  The  symbolic  form  for  the 
special  case  of  the  disk  located  on  and  aligned  with  the  projectile  axis  of  symmetry  is  also 
included  in  Appendix  B. 


The  four  roots  of  the  characteristic  equation  shown  in  equation  69  are  given  in  equation  70: 


^(^A  +  E-iF±  ■7(A  -  Ef  +  4CD  -  +  2i{AF  -  2BD  -  EF) 

if  A  +  E  +  /F  ±  V(A  -  e)"  +  ACD  -F^  +  2i{EF  +  2BD  -  AF) 


(70) 
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These  results  are  identical  to  conventional  rigid  projectile  analysis.  Consequently,  rotating 
internal  part  projectile  analysis  can  be  approached  in  essentially  the  same  manner  that  rigid 
projectiles  are  analyzed. 


5.  Example  Results 


In  order  to  examine  the  changes  an  internal  rotating  disk  induces  on  the  dynamic  behavior  of  a 
projectile,  the  following  analysis  documents  how  the  fast  and  slow  epicyclic  modes  change  for 
various  rotating  disk  arrangements.  Results  are  shown  for  a  typical  155-mm  spin-stabilized 
artillery  shell  having  a  nominal  weight  of  94.88  Ibf  with  a  reference  area  of  5  =  0.20  ft  and  a 
reference  diameter  of  D  =  0.51  ft.  The  nominal  stationline,  buttline,  and  waterline  distances  to  the 
projectile’s  center  of  gravity  are  =  1 .06  ft,  =  0.00  ft,  and  =  0.00  ft.  The  standard 

aerodynamic  center  of  pressure  and  the  aerodynamic  center  of  Magnus  are  located  at  the  following 
distances  along  the  stationline,  buttline,  and  waterline:  1.76  ft,  =  0.00  ft, 

=  0.00  ft,  and  =  1.76  ft,  =  0.00  ft,  =  0.00  ft.  The  moments  of  inertia  about  the  body- 
axis  are:  Iexx  =  0.11  slug  -  ft^,  and  Ieyy  =  1  -40  slug  -  ft^,  and  Iezz  =  1  -40  slug  -  ft^.  The  nominal 
disk  is  0.33  ft  in  diameter,  0.06  ft  thick,  and  weighs  9.5  Ibf.  It  is  nominally  located  on  the 
projectile  center  of  gravity  and  has  the  following  inertia  properties:  Idxx  =  0.0041  slug  -  ft^,  lory  = 
0.0021  slug  -  ft^,  and  Idzz  =  0.0021  slug  -  ft^.  Euler  angles  of  the  projectile  are  0.00  rad,  0  = 
0.22  rad,  and  \i/=  0.00  rad.  The  projectile  has  a  forward  velocity  of  m  =  2710  ft/s  and  a  spin  rate 
of  p  =  1674.10  rad/s.  The  pitch  and  yaw  rates  and  the  side  velocities  are  all  equal  to  zero.  The 

disk  spin  rate  and  angle  are  varied. 

In  the  analysis  to  follow,  the  ratio  of  the  disk  mass  to  projectile  mass  is  dubbed  the  mass  ratio, 
Mr,  the  ratio  of  the  disk  spin  rate  to  the  projectile  spin  rate  is  called  the  spin  ratio,  Sr,  and  the 
angle  that  the  disk  spin  axis  is  precessed  from  the  projectile  axis  is  given  as  the  disk  angle,  (///). 
The  disk  angle  is  defined  in  the  inset  of  Figure  4  for  a  disk  located  off  the  projectile  axis  of 
symmetry  in  the  buttline  direction  ( J ^  axis).  A  disk  located  on  the  projectile  axis  of  symmetry 
can  be  visualized  by  setting  the  buttline  offset  distance  equal  to  zero. 

The  effect  on  the  epicyclic  dynamics  of  the  system  of  varying  the  disk  angle  for  three  different 
mass  ratios  with  a  spin  ratio  of  5  is  demonstrated  in  the  root- locus  plots  shown  in  Figures  4-6. 
The  disk  angle  is  varied  from  0°  (represented  by  the  diamond)  to  360°.  The  circle  represents  the 
eigenvalues  for  a  disk  angle  of  90°,  the  square  for  a  disk  angle  of  180°,  and  the  bold  X  denotes 
the  location  of  the  eigenvalues  determined  from  a  similar  rigid  artillery  round  without  a  rotating 
internal  part.  These  conventions  are  used  for  all  root-locus  plots.  In  Figures  4-8,  the  epicyclic 
modes  travel  from  the  diamond  to  the  square  as  the  disk  angle  is  increased  from  0°  to  180°,  and 
march  back  up  the  same  path  towards  the  diamond  as  the  disk  angle  is  further  increased  from 
180°  to  360°.  In  Figure  4  it  can  be  seen  that  a  disk  of  mass  ratio  I/IOO  slightly  affects  the 
projectile  fast  and  slow  modes  at  angles  other  than  90°.  Comparison  of  the  three  root- locus  plots 
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Figure  8.  Root-locus  for  spin  ratio  =  1  and  mass  ratio  =  1/10  (0  =  0°,  360°  disk  angle; 
o  =  90°  disk  angle;  □  =  180°  disk  angle;  X  =  rigid  projectile). 

reveals  that  increasing  the  mass  ratio  increases  the  effect  that  the  spinning  disk  has  on  the 
epicyclic  dynamics.  As  shown  in  Figure  6,  when  the  mass  ratio  is  sufficiently  large,  the  fast 
mode  can  become  unstable  for  large  disk  angles  greater  than  90°.  In  this  case,  spin  stability  of 
the  complete  round  is  adversely  impacted  by  the  component  of  disk  angular  velocity  in  the 
opposite  direction  of  projectile  spin.  Also  note  that  the  mass  ratio  has  very  little  effect  on  the 
stability  of  the  system  for  a  disk  angle  of  90°.  For  disk  angles  >180°,  the  eigenvalues  march 
back  along  the  same  path  toward  0°.  This  is  expected  because  for  a  disk  located  on  the  projectile 
axis  of  symmetry  the  disk  angles  from  0°  to  180°  relative  to  the  projectile  are  equivalent  to  the 
angles  from  180°  to  360°  with  the  direction  of  spin  reversed. 

Figures  7-9  are  the  root- locus  plots  obtained  by  varying  the  disk  angle  for  three  different  spin 
ratios  with  a  disk  to  projectile  mass  ratio  of  1/10.  By  comparing  Figures  7-9,  it  is  shown  that 
decreasing  spin  ratio  diminishes  the  effect  that  a  spinning  disk  of  a  given  mass  has  on  the 
epicyclic  modes.  The  opposite  holds  true  as  well.  Increasing  the  spin  ratio  increases  the 
dynamic  effects  and  is  capable  of  driving  the  system  unstable  for  certain  disk  angles. 

Comparison  of  Figures  6  and  9  shows  that  the  same  dynamic  effects  can  be  achieved  with  either 
mass  ratio  or  spin  ratio.  Figures  4-9  also  demonstrate  that  a  disk  of  any  mass  ratio,  spinning  in 
the  same  direction  as  the  projectile  (i.e.,  disk  angle  <90°),  tends  to  stabilize  the  fast  mode  of  the 
system  while  destabilizing  the  slow  mode.  Disk  angles  of  >90°  up  to  180°  have  the  opposite 
effect.  However,  for  large  spin  ratios  and  large  mass  ratios,  the  movement  of  the  modes  is  much 
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Figure  9.  Root-locus  for  spin  ratio  =  10  and  mass  ratio  =  1/10  (0  =  0°,  360°  disk  angle; 
o  =  90°  disk  angle;  □  =  180°  disk  angle;  X  =  rigid  projectile). 

greater  per  increase  in  disk  angle.  While  not  explicitly  shown  previously,  stationline  position  of 
the  disk  has  no  effect  on  the  fast  and  slow  modes  when  the  disk  mass  center  is  located  on  the 
projectile  axis  of  symmetry.  Moreover,  in  this  case,  if  the  product  MrSr  is  held  constant,  the  fast 
and  slow  modes  do  not  change  with  disk  location  along  the  projectile  axis  of  symmetry. 

Figures  9  and  10  examine  the  fast  and  slow  modes  for  spinning  disks  that  are  located  off  the 
projectile  axis  of  symmetry  in  the  buttline  direction  and  the  waterline  direction.  It  is  important  to 
make  this  distinction  because  rotating  the  disk  about  the  Kg  axis  from  the  projectile  axis  of 
symmetry  produces  two  unique  angles  relative  to  the  projectile  in  each  direction.  Analysis  of  the 
fast  and  slow  modes  for  a  system  with  the  disk  located  off  the  projectile  axis  of  symmetry  in  the 
buttline  direction  and  varying  its  position  in  the  direction  of  the  stationline  produce  the  same 
mode  path  as  that  for  a  disk  located  on  the  projectile  axis  of  symmetry  and  varying  its  stationline 
location.  The  path  is  the  same,  however,  the  locations  of  the  modes  along  the  path  do  vary 
slightly  with  stationline  position.  This  slight  variation  is  undetectable  on  the  scale  used  for 
graphs  shown  in  this  report.  Therefore,  Figure  9  can  be  referred  to  for  a  root- locus  of  the 
previously  mentioned  system  for  a  disk  located  anywhere  along  the  stationline  and  anywhere  off 
the  axis  in  the  buttline  direction.  The  stationline  location  does,  however,  change  the  fast  and 
slow  modes  of  the  system  for  a  disk  that  is  located  off-axis  in  the  waterline  direction  as  shown  in 
Figure  10. 
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Figure  10.  Root-locus  for  a  disk  located  3/4  in  off-axis  of  symmetry  in  the  waterline  direction  (spin  ratio  =  10, 

mass  ratio  =  1/10)  (0  =  0°,  360°  disk  angle;  o  =  90°  disk  angle;  □  =  180°  disk  angle;  X  =  rigid  projectile). 

In  Figure  10,  the  fast  and  slow  modes  are  plotted  in  increments  of  disk  angle,  and  to  add  clarity 
to  the  graphs  only  the  eigenvalues  containing  positive  imaginary  parts  are  shown.  The  root-locus 
for  the  stationline  of  the  disk  center  of  gravity  that  coincides  with  the  stationline  of  the  projectile 
center  of  gravity  is  the  same  as  that  for  disks  either  located  on  the  axis  of  symmetry  or  off-axis  in 
the  buttline  direction.  The  plots  are  only  shown  for  stationline  positions  varied  from  the  aft  of 
the  projectile  to  its  center  of  gravity.  Disks  that  are  located  forward  of  the  projectile  center  of 
gravity,  at  distances  from  the  center  of  gravity  equal  to  those  that  are  aft  of  it,  produce  the  same 
modes  for  the  respective  distances.  Notice  that  the  fast  modes  for  SL  =  0.0627  ft  and 
SL  =  0.5627  ft  become  slow  modes  at  a  disk  angle  of  180°  (designated  by  the  square)  and  the 
slow  modes  change  to  fast  modes.  It  is  also  interesting  that  for  these  two  curves  the  directions  of 
the  root-locus  path  change  for  disk  angles  slightly  <180°  and  slightly  >90°  (designated  by  the 
circle).  The  previously  mentioned  studies  show  that  for  a  disk  located  at  the  projectile  center  of 
gravity  stationline,  the  epicyclic  dynamics  of  the  system  are  not  affected  by  the  off-axis  disk 
location.  They  also  show  that  a  spinning  disk  mounted  at  an  angle  of  90°  with  a  stationline 
position  equal  to  the  stationline  of  the  center  of  gravity  of  the  projectile  is  the  same  system  as  a 
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rigid  projectile  without  a  disk.  The  slow  mode  of  the  system  with  a  disk  located  off-axis  in  the 
buttline  direction  is  predominantly  on  the  left  side  of  the  s-plane,  whereas  for  a  disk  located  off- 
axis  in  the  direction  of  the  waterline  the  slow  mode  is  divided  between  the  left  and  right  sides  of 
the  s-plane. 

Figures  1 1  and  12  show  the  root- locus  for  various  angles  of  a  disk  located  off-axis  in  various 
directions  lying  between  the  waterline  and  buttline.  The  fast  mode  shown  in  Figure  1 1  and  the 
slow  mode  in  Figure  12  are  for  a  disk  located  1/4  in  off  the  projectile  axis  of  symmetry  with  a 
spin  ratio  of  10  and  mass  ratio  of  1/10  mounted  at  0.0627  ft  along  the  stationline.  The  lines 
running  vertically  represent  the  epicyclic  mode  locations  for  different  directions  of  offset  defined 
in  the  figure  legend.  The  center  of  gravity  location  is  the  same  regardless  of  the  disk  angle.  The 
vertical  line  running  roughly  down  the  center  of  the  plot  represents  the  epicyclic  modes  for  a  disk 
located  solely  in  the  positive  buttline  direction.  Each  horizontal  line  in  the  graph  represents  the 
epicyclic  mode  locations  for  a  specific  disk  angle  designated  by  the  values  shown.  The 
horizontal  lines  located  on  the  left-hand  side  of  the  center  line  represent  angles  from  0°  to  180°, 
and  those  on  the  right-hand  side  represent  angles  from  180°  to  360°.  If  the  fast  mode  of  a  disk 
spinning  at  an  angle  of  120°  and  located  45°  from  the  positive  buttline  direction  were  desired  to 
be  known,  one  would  follow  the  dash-dot  line  until  it  crossed  the  bottom  half  of  the  120°  curve 
and  arrive  at  a  value  of  -0.53E  ^-i-0.017/.  These  plots  definitively  show  that  the  epicyclic  modes 
of  a  projectile  containing  a  rotating  internal  part  are  significantly  affected  by  the  location  and 
angle  of  the  rotating  internal  part.  Eigure  13  shows  a  system  with  a  disk  mounted  at  0.0627  ft 
along  the  stationline  and  at  various  distances  off-axis  in  a  direction  equidistant  from  the  buttline 
and  waterline.  Increasing  the  distance  that  the  disk  is  located  off-axis  increases  the  effects  that 
the  spinning  disk  has  on  the  dynamics  of  the  system.  Greater  distances  off-axis  cause  the  fast 
and  slow  modes  to  swap  at  disk  angles  near  180°,  whereas  smaller  distances,  such  as  1/4  in, 
shown  in  Eigures  1 1  and  12,  do  not.  Eigure  14  shows  the  root- locus  for  a  disk  located  3/4  in  off 
the  projectile  axis  of  symmetry  in  a  direction  equidistant  from  the  waterline  and  buttline 
directions  for  various  spin  ratios.  Eike  disks  located  on  the  axis  of  symmetry,  increasing  the  spin 
ratio  increases  the  range  of  the  epicyclic  modes.  Eor  disk  angles  between  0°  and  180°  the  fast 
mode  is  predominantly  on  the  left-hand  side  of  the  imaginary  axis,  and  the  slow  mode  is 
predominantly  on  the  right-hand  side.  The  opposite  is  true  for  disk  angles  between 
180°  and  360°.  Increasing  the  spin  ratio  moves  the  modes  farther  from  the  imaginary  axis  for 
disk  angles  other  than  0°  and  360°  and  expands  their  range  along  the  imaginary  axis. 
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Figure  13.  Root-locus  for  a  disk  located  at  various  distances  off-axis  of  symmetry 
in  a  direction  equidistant  from  the  waterline  and  buttline  directions  (spin 
ratio  =  10;  mass  ratio  =  1/10)  (0  =  0°,  360°  disk  angle;  o  =  90°  disk 
angle;  □  =  180°  disk  angle;  X  =  rigid  projectile). 


Figure  14.  Root-locus  for  a  disk  located  3/4  in  off-axis  of  symmetry  in  a  direction 

equidistant  from  the  waterline  and  buttline  directions  for  various  spin  ratios 
(mass  ratio  =  1/10)  (0  =  0°,  360°  disk  angle;  o  =  90°  disk  angle;  □  =  180° 
disk  angle;  X  =  rigid  projectile). 
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6.  Conclusions 


The  equations  of  motion  for  a  projectile  containing  an  axisymmetric  rotating  internal  disk  that 
spins  at  a  constant  rate  have  been  developed.  The  model  allows  for  the  disk  to  be  located  off  the 
axis  of  symmetry  of  the  projectile  and  oriented  at  an  arbitrary  angle  relative  to  the  projectile  axis 
of  symmetry.  Projectile  linear  theory  has  been  modified  to  accommodate  projectile 
configurations  that  contain  an  internal  rotating  disk.  The  addition  of  an  internal  rotating  disk 
alters  several  of  the  coefficients  in  the  epicyclic  dynamic  equations  leading  to  modified  fast  and 
slow  epicyclic  modes. Using  modified  projectile  linear  theory,  the  effect  of  disk  mass,  speed,  and 
location  is  systematically  studied.  If  the  mass  ratio  times  the  spin  ratio  is  held  constant  the  same 
epicyclic  modes  are  produced  for  a  disk  located  at  any  position  along  the  projectile  axis  of 
symmetry.  For  a  specified  mass  and  spin  ratio  the  stationline  position  of  the  disk  has  no  effect 
on  the  dynamics  of  the  system  for  a  disk  located  either  on  the  projectile  axis  of  symmetry. 
However,  the  stationline  position  does  affect  the  dynamics  if  the  disk  is  off-axis.  In  these 
instances,  the  farther  the  disk  is  located  away  from  the  center  of  gravity  of  the  projectile  the  more 
influence  its  mass,  spin  rate,  and  angle  relative  to  the  projectile  axis  of  symmetry  have  on  the 
epicyclic  dynamics  of  the  system. 
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Appendix  A.  Coefficients  of  the  Epicyclic  Dynamic  Equations  in  C-Form 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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The  coefficients  of  the  characteristic  equation  in  C-programming  language  form  are  given  in 
equations  A1-A6.  The  equations  were  derived  using  Wolfram  Research  Inc.’s  Mathematica  and 
the  following  variable  substitutions  need  to  be  made  in  order  for  their  appearance  to  be 
consistent  with  the  text. 


dx,dy ,  dz  <=>  x ^  y ^  z 
phi,  theta,  psi  <::>(j),6,y/ 

CYPA  Cypj, 

CLP  Cpp 

CXI  <s> 
omega  Q. 


Variable  Substitutions 

dax,  day,  daz  ^PA-^y pa-’^pa 

md  nip, 

CNA  « 

CDD  Coo 

CYB1,CZB1  <z>Cyp„CzBi 
Cos(disk),  Sin(disk), 

0  ,ny  ,n2 


dmx,  dmy,  dmz 

^  ^PM  ’  y  PM  ’  ^PM 

mp  nip 
CMQ  ^ 

CXO,  CYO,  CZO 

^  Cxo’  ^yo’  ^zo 
rho-^  p 


pxx  ’  Cyy  ’  Czz  ’  Cxv  ’  Cxz  ’  Cyz 


JPXX,  JPYY,  JPZZ,  JPXY,  JPXZ,  JPYZ  « I 
JDXX,  JDYY,  JDZZ,  JDXY,  JDXZ,  JDYZ  ,  I ^yy ,  Ip^^z ,  Cxy  .  Cxz .  Cyz 


Al)  A=  -(CNA*D*rho*S)/(2.*mc) 


A2)  B  =  (D*(-((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*(-(CYPA*D*dy*md*pt*rho*S*Cos(phi))/(4.*(md  +  mp))  +  & 
(CNA*dz*md*rho*S*V*Cos(phi))/(2.*(md  +  mp))  +  & 
(CYPA*D*dz*md*pt*rho*S*Sin(phi))/(4.*(md  +  mp))  +  & 
(CNA*dy*md*rho*S*V*Sin(phi))/(2.*(md  +  mp))))  +  (JPXX  +  ((Power(dy,2)  + 

& 

Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CYPA*D*dmx*pt*rho*S)/4.  +  & 
(CYPA*D*dx*md*pt*rho*S*Power(Cos(phi),2))/(4.*(md  +  mp))  +  & 
(CYPA*D*dx*md*pt*rho*S*Power(Sin(phi),2))/(4.*(md  +  & 
inp))))/(-Power(Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2) 
+  & 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  & 

JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
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mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp)))  +  & 
((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((-((Cos(phi)*(JPXY  - 
(dx*dy*md*mp)/(md  & 

+  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md 
+  & 

mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md 
+  & 

mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  + 
& 

mp))*Sin(phi))*(-(CYPA*D*dy*md*pt*rho*S*Cos(phi))/(4.*(md  +  mp))  +  & 
(CNA*dz*md*rho*S*V*Cos(phi))/(2.*(md  +  mp))  +  & 
(CYPA*D*dz*md*pt*rho*S*Sin(phi))/(4.*(md  +  mp))  +  & 
(CNA*dy*md*rho*S*V*Sin(phi))/(2.*(md  +  mp))))  +  (JPXX  +  ((Power(dy,2)  + 

& 

Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CYPA*D*dmx*pt*rho*S)/4.  +  & 
(CYPA*D*dx*md*pt*rho*S*Power(Cos(phi),2))/(4.*(md  +  mp))  +  & 
(CYPA*D*dx*md*pt*rho*S*Power(Sin(phi),2))/(4.*(md  +  mp)))))  +  & 
(-((-(CYPA*D*dy*md*pt*rho*S*Cos(phi))/(4.*(md  +  mp))  +  & 
(CNA*dz*md*rho*S*V*Cos(phi))/(2.*(md  +  mp))  +  & 
(CYPA*D*dz*md*pt*rho*S*Sin(phi))/(4.*(md  +  mp))  +  & 
(CNA*dy*md*rho*S*V*Sin(phi))/(2.*(md  +  mp)))*((JPXZ  - 
(dx*dz*md*mp)/(md  +  & 

mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
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JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*((CNA*dax*rho*S*V)/2.  -  & 
(CNA*dx*md*rho*S*V*Power(Cos(phi),2))/(2.*(md  +  mp))  -  & 
(CNA*dx*md*rho*S*V*Power(Sin(phi),2))/(2.*(md  +  mp))))*(- 
Power(Cos(phi)*(JPXY  & 

-  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))))/((-Power(Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2) 
+  & 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  & 

JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))*(-Power(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp),2)  +  (-Power(Cos(phi)*(JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
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JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))*(-Power((JPXZ  -  & 
(dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  + 
(JDXX-  & 

JDYY)*Cos(disk)*Sin(disk))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md 
+  mp)  & 

+  JPYY*(md  +  mp))*Power(Cos(phi),2)  -  2*dy*dz*md*mp*Cos(phi)*Sin(phi)  + 
& 

(JPYY*md  +  JPYY*mp  +  Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  + 
JDYY*(md+  & 

mp)*Power(Cos(disk),2)  +  JDXX*(md  + 
mp) *Power(Sin(disk) ,2)) *Power(Sin(phi),2)  & 

+  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))))))/V 


A3)  C  =  (D*(-((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*(-(CYPA*D*dz*md*pt*rho*S*Cos(phi))/(4.*(md  +  mp))  -  & 
(CNA*dy*md*rho*S*V*Cos(phi))/(2.*(md  +  mp))  -  & 
(CYPA*D*dy*md*pt*rho*S*Sin(phi))/(4.*(md  +  mp))  +  & 
(CNA*dz*md*rho*S*V*Sin(phi))/(2.*(md  +  mp))))  +  (JPXX  +  ((Power(dy,2)  + 

& 

Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CNA*dax*rho*S*V)/2.  +  & 
(CNA*dx*md*rho*S*V*Power(Cos(phi),2))/(2.*(md  +  mp))  +  & 
(CNA*dx*md*rho*S*V*Power(Sin(phi),2))/(2.*(md  +  mp))))/(- 
Power(Cos(phi)*(JPXY  & 

-  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp)))  +  & 
((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
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JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((-((Cos(phi)*(JPXY  - 
(dx*dy*md*mp)/(md  & 

+  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md 
+  & 

mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md 
+  & 

mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  + 
& 

mp))*Sin(phi))*(-(CYPA*D*dz*md*pt*rho*S*Cos(phi))/(4.*(md  +  mp))  -  & 
(CNA*dy*md*rho*S*V*Cos(phi))/(2.*(md  +  mp))  -  & 
(CYPA*D*dy*md*pt*rho*S*Sin(phi))/(4.*(md  +  mp))  +  & 
(CNA*dz*md*rho*S*V*Sin(phi))/(2.*(md  +  mp))))  +  (JPXX  +  ((Power(dy,2)  + 

& 

Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CNA*dax*rho*S*V)/2.  +  & 
(CNA*dx*md*rho*S*V*Power(Cos(phi),2))/(2.*(md  +  mp))  +  & 
(CNA*dx*md*rho*S*V*Power(Sin(phi),2))/(2.*(md  +  mp)))))  +  & 
(-((-(CYPA*D*dz*md*pt*rho*S*Cos(phi))/(4.*(md  +  mp))  -  & 
(CNA*dy*md*rho*S*V*Cos(phi))/(2.*(md  +  mp))  -  & 
(CYPA*D*dy*md*pt*rho*S*Sin(phi))/(4.*(md  +  mp))  +  & 
(CNA*dz*md*rho*S*V*Sin(phi))/(2.*(md  +  mp)))*((JPXZ  - 
(dx*dz*md*mp)/(md  +  & 

mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
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JDYY*Power(Sin(disk),2))*(-(CYPA*D*dmx*pt*rho*S)/4.  +  & 
(CYPA*D*dx*md*pt*rho*S*Power(Cos(phi),2))/(4.*(md  +  mp))  +  & 
(CYPA*D*dx*md*pt*rho*S*Power(Sin(phi),2))/(4.*(md  +  & 
mp))))*(-Power(Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2) 
+  & 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  & 

JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))))/((-Power(Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2) 
+  & 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  & 

JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))*(-Power(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp),2)  +  (-Power(Cos(phi)*(JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
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mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))*(-Power((JPXZ  -  & 
(dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  + 
(JDXX-  & 

JDYY)*Cos(disk)*Sin(disk))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md 
+  mp)  & 

+  JPYY*(md  +  mp))*Power(Cos(phi),2)  -  2*dy*dz*md*mp*Cos(phi)*Sin(phi)  + 
& 

(JPYY*md  +  JPYY*mp  +  Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  + 
JDYY*(md+  & 

mp)*Power(Cos(disk),2)  +  JDXX*(md  + 
mp)  *Power(Sin(disk)  ,2))*Power(Sin(phi)  ,2)  & 

+  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))))))/V 


A4)  D  =  D 


A5)  E  =  (D*(-((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 

JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*(JDXX*omega*Power(Cos(disk),2)*Sin(disk)*Sin(phi)  +  & 
JDXX*omega*Power(Sin(disk),3)*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CMQ*Power(D,2)*rho*S*V)/4.  -  & 
2*JPYZ*pt*Power(Cos(phi),2)  +  (2*dy*dz*md*mp*pt*Power(Cos(phi),2))/(md  + 
mp)  & 

+  2*JDYY*pt*Cos(phi)*Sin(phi)  +  & 

(2*Power(dy,2)*md*mp*pt*Cos(phi)*Sin(phi))/(md  +  mp)  -  & 
(2*Power(dz,2)*md*mp*pt*Cos(phi)*Sin(phi))/(md  +  mp)  -  & 
2*JDYY*pt*Power(Cos(disk),2)*Cos(phi)*Sin(phi)  -  & 
2*JDXX*pt*Cos(phi)*Power(Sin(disk),2)*Sin(phi)  + 
2*JPYZ*pt*Power(Sin(phi),2)  & 

-  (2*dy*dz*md*mp*pt*Power(Sin(phi),2))/(md  +  mp)))/(-Power(Cos(phi)*(JPXY 

-  & 

(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*  Cos  (disk)  *Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 
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JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp)))  +  & 
((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((-((Cos(phi)*(JPXY  - 
(dx*dy*md*mp)/(md  & 

+  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md 
+  & 

mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md 
+  & 

mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  + 
& 

mp))*Sin(phi))*(JDXX*omega*Power(Cos(disk),2)*Sin(disk)*Sin(phi)  +  & 
JDXX*omega*Power(Sin(disk),3)*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CMQ*Power(D,2)*rho*S*V)/4.  -  & 
2*JPYZ*pt*Power(Cos(phi),2)  +  (2*dy*dz*md*mp*pt*Power(Cos(phi),2))/(md  + 
mp)  & 

+  2*JDYY*pt*Cos(phi)*Sin(phi)  +  & 

(2*Power(dy,2)*md*mp*pt*Cos(phi)*Sin(phi))/(md  +  mp)  -  & 
(2*Power(dz,2)*md*mp*pt*Cos(phi)*Sin(phi))/(md  +  mp)  -  & 
2*JDYY*pt*Power(Cos(disk),2)*Cos(phi)*Sin(phi)  -  & 
2*JDXX*pt*Cos(phi)*Power(Sin(disk),2)*Sin(phi)  + 
2*JPYZ*pt*Power(Sin(phi),2)  & 
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-  (2*dy*dz*md*mp*pt*Power(Sin(phi),2))/(md  +  mp))))  +  & 
(-((JDXX*omega*Power(Cos(disk),2)*Sin(disk)*Sin(phi)  +  & 
JDXX*omega*Power(Sin(disk),3)*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  & 
mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(JDYY*pt*Power(Cos(phi),2))  -  & 
JPXX*pt*Power(Cos(phi),2)  - 

(2*Power(dy,2)*md*mp*pt*Power(Cos(phi),2))/(md  +  & 
mp)  -  JDXX*pt*Power(Cos(disk),2)*Power(Cos(phi),2)  +  & 
JDYY*pt*Power(Cos(disk),2)*Power(Cos(phi),2)  -  & 
JDXX*omega*Power(Cos(disk),3)*Power(Cos(phi),2)  +  & 
JDXX*pt*Power(Cos(phi),2)*Power(Sin(disk),2)  -  & 
JDYY*pt*Power(Cos(phi),2)*Power(Sin(disk),2)  -  & 
JDXX*omega*Cos(disk)*Power(Cos(phi),2)*Power(Sin(disk),2)  -  & 
4*JPYZ*pt*Cos(phi)*Sin(phi)  +  (4*dy*dz*md*mp*pt*Cos(phi)*Sin(phi))/(md  + 
mp)  & 

+  JDYY*pt*Power(Sin(phi),2)  -  JPXX*pt*Power(Sin(phi),2)  -  & 
(2*Power(dz,2)*md*mp*pt*Power(Sin(phi),2))/(md  +  mp)  -  & 
JDXX*pt*Power(Cos(disk),2)*Power(Sin(phi),2)  -  & 
JDYY*pt*Power(Cos(disk),2)*Power(Sin(phi),2)  -  & 
JDXX*omega*Power(Cos(disk),3)*Power(Sin(phi),2)  -  & 
JDXX*pt*Power(Sin(disk),2)*Power(Sin(phi),2)  -  & 
JDYY*pt*Power(Sin(disk),2)*Power(Sin(phi),2)  -  & 
JDXX*omega*Cos(disk)*Power(Sin(disk),2)*Power(Sin(phi),2)))*(- 
Power(Cos(phi)*(  & 

JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (- 
JPXZ  & 

+  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))))/((-Power(Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2) 
+  & 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  & 

JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
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& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))*(-Power(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp),2)  +  (-Power(Cos(phi)*(JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))*(-Power((JPXZ  -  & 
(dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  + 
(JDXX  -  & 

JDYY)*Cos(disk)*Sin(disk))*Sin(phi),2)  -i-  ((JPXX  -i-  ((Power(dy,2)  -i-  & 
Power(dz,2))*md*mp)/(md  -i-  mp)  -i-  JDXX*Power(Cos(disk),2)  -i-  & 
JDYY*Power(Sin(disk),2))*(((Power(dx,2)  -i-  Power(dy,2))*md*mp  -i-  JDYY*(md 
-I-  mp)  & 

-I-  JPYY*(md  -I-  mp))*Power(Cos(phi),2)  -  2*dy*dz*md*mp*Cos(phi)*Sin(phi)  -i- 
& 

(JPYY*md  -I-  JPYY*mp  -i-  Power(dx,2)*md*mp  -i-  Power(dz,2)*md*mp  -i- 
JDYY*(md-i-  & 

mp)*Power(Cos(disk),2)  -i-  JDXX*(md  -i- 
mp)  *Power(Sin(disk)  ,2))*Power(Sin(phi)  ,2)  & 

-I-  JPYZ*(md  -I-  mp)*Sin(2*phi)))/(md  -i-  mp))))))/V 


A6)  F  =  -((D*(-((-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  -i-  mp)  -i-  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))  -i-  (-JPXZ  -i-  (dx*dz*md*mp)/(md  -i-  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  -i-  mp))*Cos(phi)  -i-  (JPXY  -  & 
(dx*dy*md*mp)/(md  -i-  mp)  -i-  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  -i- 
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& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp))*(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md 
+  & 

mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  + 
& 

mp))*Sin(phi))*(JDXX*omega*Power(Cos(disk),2)*Sin(disk)*Sin(phi)  +  & 
JDXX*omega*Power(Sin(disk),3)*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(CMQ*Power(D,2)*rho*S*V)/4.  -  & 
2*JPYZ*pt*Power(Cos(phi),2)  +  (2*dy*dz*md*mp*pt*Power(Cos(phi),2))/(md  + 
mp)  & 

+  2*JDYY*pt*Cos(phi)*Sin(phi)  +  & 

(2*Power(dy,2)*md*mp*pt*Cos(phi)*Sin(phi))/(md  +  mp)  -  & 
(2*Power(dz,2)*md*mp*pt*Cos(phi)*Sin(phi))/(md  +  mp)  -  & 
2*JDYY*pt*Power(Cos(disk),2)*Cos(phi)*Sin(phi)  -  & 
2*JDXX*pt*Cos(phi)*Power(Sin(disk),2)*Sin(phi)  + 
2*JPYZ*pt*Power(Sin(phi),2)  & 

-  (2*dy*dz*md*mp*pt*Power(Sin(phi),2))/(md  +  mp))))  +  & 
(-((JDXX*omega*Power(Cos(disk),2)*Sin(disk)*Sin(phi)  +  & 
JDXX*omega*Power(Sin(disk),3)*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  & 
mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  & 
JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  +  (JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(-(JDYY*pt*Power(Cos(phi),2))  -  & 
JPXX*pt*Power(Cos(phi),2)  - 

(2*Power(dy,2)*md*mp*pt*Power(Cos(phi),2))/(md  +  & 
mp)  -  JDXX*pt*Power(Cos(disk),2)*Power(Cos(phi),2)  +  & 
JDYY*pt*Power(Cos(disk),2)*Power(Cos(phi),2)  -  & 
JDXX*omega*Power(Cos(disk),3)*Power(Cos(phi),2)  +  & 
JDXX*pt*Power(Cos(phi),2)*Power(Sin(disk),2)  -  & 
JDYY*pt*Power(Cos(phi),2)*Power(Sin(disk),2)  -  & 
JDXX*omega*Cos(disk)*Power(Cos(phi),2)*Power(Sin(disk),2)  -  & 
4*JPYZ*pt*Cos(phi)*Sin(phi)  +  (4*dy*dz*md*mp*pt*Cos(phi)*Sin(phi))/(md  + 
mp)  & 

+  JDYY*pt*Power(Sin(phi),2)  -  JPXX*pt*Power(Sin(phi),2)  -  & 
(2*Power(dz,2)*md*mp*pt*Power(Sin(phi),2))/(md  +  mp)  -  & 
JDXX*pt*Power(Cos(disk),2)*Power(Sin(phi),2)  -  & 
JDYY*pt*Power(Cos(disk),2)*Power(Sin(phi),2)  -  & 
JDXX*omega*Power(Cos(disk),3)*Power(Sin(phi),2)  -  & 
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JDXX*pt*Power(Sin(disk),2)*Power(Sin(phi),2)  -  & 
JDYY*pt*Power(Sin(disk),2)*Power(Sin(phi),2)  -  & 
JDXX*omega*Cos(disk)*Power(Sin(disk),2)*Power(Sin(phi),2)))*(- 
Power(Cos(phi)*(  & 

JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (- 
JPXZ  & 

+  (dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  & 
mp))))/(V*(-Power(-((Cos(phi)*(JPXY  -  (dx*dy*md*mp)/(md  +  mp)  +  (JDXX  - 
& 

JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  +  (dx*dz*md*mp)/(md  +  & 
mp))*Sin(phi))*((JPXZ  -  (dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))*Sin(phi)))  + 
& 

((JPXX  +  ((Power(dy,2)  +  Power(dz,2))*md*mp)/(md  +  mp)  +  & 
JDXX*Power(Cos(disk),2)  +  JDYY*Power(Sin(disk),2))*((-(dy*dz*md*mp)  +  & 
JPYZ*(md  +  mp))*Power(Cos(phi),2)  +  Cos(phi)*(-(JDYY*md)  -  JDYY*mp  -  & 
Power(dy,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))*Sin(phi)  -  (-(dy*dz*md*mp)  +  JPYZ*(md 
+  & 

mp))*Power(Sin(phi),2)))/(md  +  mp),2)  +  (-Power(Cos(phi)*(JPXY  -  & 
(dx*dy*md*mp)/(md  +  mp)  +  (JDXX  -  JDYY)*Cos(disk)*Sin(disk))  +  (-JPXZ  + 
& 

(dx*dz*md*mp)/(md  +  mp))*Sin(phi),2)  +  ((JPXX  +  ((Power(dy,2)  +  & 
Power(dz,2))*md*mp)/(md  +  mp)  +  JDXX*Power(Cos(disk),2)  +  & 
JDYY*Power(Sin(disk),2))*(Power(Cos(phi),2)*(JPYY*md  +  JPYY*mp  +  & 
Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  +  JDYY*(md  + 
mp)*Power(Cos(disk),2)  +  & 

JDXX*(md  +  mp)*Power(Sin(disk),2))  +  Sin(phi)*(2*dy*dz*md*mp*Cos(phi)  + 
& 

((Power(dx,2)  +  Power(dy,2))*md*mp  +  JDYY*(md  +  mp)  +  JPYY*(md  +  & 
mp))*Sin(phi))  -  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))*(-Power((JPXZ  -  & 
(dx*dz*md*mp)/(md  +  mp))*Cos(phi)  +  (JPXY  -  (dx*dy*md*mp)/(md  +  mp)  + 
(JDXX  -  & 

JDYY)*Cos(disk)*Sin(disk))*Sin(phi),2)  -i-  ((JPXX  -i-  ((Power(dy,2)  -i-  & 
Power(dz,2))*md*mp)/(md  -i-  mp)  -i-  JDXX*Power(Cos(disk),2)  -i-  & 
JDYY*Power(Sin(disk),2))*(((Power(dx,2)  -i-  Power(dy,2))*md*mp  -i-  JDYY*(md 
-I-  mp)  & 
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+ 


JPYY*(md  +  mp))*Power(Cos(phi),2)  -  2*dy*dz*md*mp*Cos(phi)*Sin(phi)  + 

(JPYY*md  +  JPYY*mp  +  Power(dx,2)*md*mp  +  Power(dz,2)*md*mp  + 
JDYY*(md+  & 

mp)*Power(Cos(disk),2)  +  JDXX*(md  + 
mp)  *Power(Sin(disk)  ,2))*Power(Sin(phi)  ,2)  & 

+  JPYZ*(md  +  mp)*Sin(2*phi)))/(md  +  mp))))) 
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Appendix  B.  Coefficients  of  the  Epicyclic  Dynamic  Equations  in  Symbolic- 
Eorm  for  a  Special  Case 


Application  of  the  linear  theory  assumptions  to  the  simplified  system  of  a  projectile  containing  a 
rotating  internal  part  located  on  the  projectile  axis  of  symmetry  at  a  disk  angle,  equal  to  0° 
yields  the  coefficients  of  the  epicyclic  equations  shown  in  equations  B-1  through  B-6: 


^^  pSDC^ 
2m^ 


(B-1) 


pSD  m^D  ^PM  ^PD  nt’c  ^ypaP 

y  V  ^CYY  y  ^ 


(B-2) 


pSD  m^D  {^PA  ^PD  me  y^NA 

Jv  ^CYY 


(B-3) 


D  =  D; 

V  2m(,  If^YY  )  2 

^  ^ PXxP  ^DXX  {p 

V  )  I^yy 


(B-4) 

(B-5) 

(B-6) 


Magnus  force  is  assumed  to  be  small  in  comparison  to  other  aerodynamic  forces  and  is  dropped 
from  equations  B-1  through  B-6.  However,  due  to  the  magnitude  amplification  resulting  from 
the  cross  product  between  Magnus  force  and  its  respective  moment  arm,  the  Magnus  moment  is 
retained  in  equation  B-2.  The  equations  for  the  simplified  system  previously  presented  are 
essentially  the  same  as  those  derived  for  a  dual-spin  projectile  by  Costello  and  Peterson.^  They 
are  identical  if  the  aerodynamic  coefficients  applied  on  the  aft  section  of  a  dual-spin  projectile 
are  neglected. 


*  Costello,  M.;  Peterson,  A.  Linear  Theory  of  a  Dual  Spin  Projectile  in  Atmospheric  Flight.  Journal  of  Guidance,  Control, 


and  Dynamics  2000,  23  (4),  789-797. 
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Intentionally  lelt  blank. 
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List  of  Abbreviations  and  Symbols 


=  inertial  frame  unit  vectors 
Ijj,Jg,Kg  =  body  frame  unit  vectors 


J  J  K 


no-roll  frame  unit  vectors 


^PD’’  y  PD’’  ^PD 


Vp  , 


^PA->  y  PA  ■>  ^PA 


^PM  ^  y  PM  ^  ^PM 


u,v,w 

u,v,  w 


p,q,r 


p,q,r 


(j),0,\j/ 


x,y,z 


(Onii 


body  frame  components  of  distance  vector  from  projectile  center  of  mass 
to  disk  center  of  mass 

body  frame  components  of  distance  vector  from  projectile  base  to 
projectile  center  of  mass 

body  frame  components  of  distance  vector  from  projectile  base  to  disk 
center  of  mass 

nominal  body  frame  components  of  distance  vector  from  projectile  base  to 
projectile  center  of  mass 

no-roll  frame  components  of  distance  vector  from  projectile  center  of  mass 
to  center  of  pressure 

no-roll  frame  components  of  distance  vector  from  projectile  center  of  mass 
to  Magnus  center  of  pressure 

translational  velocity  components  of  the  two-body  system  center  of  mass 
resolved  in  the  body  frame 

translational  velocity  components  of  the  two-body  system  center  of  mass 
resolved  in  the  no-roll  frame 

roll,  pitch,  and  yaw  components  of  the  angular  velocity  vector  of  projectile 
body  expressed  in  the  body  frame 

roll,  pitch,  and  yaw  components  of  the  angular  velocity  vector  of  projectile 
body  expressed  in  the  no-roll  frame 

Euler  roll,  pitch,  and  yaw  angles 

position  vector  components  of  the  two-body  system  center  of  mass 
expressed  in  the  inertial  frame 

angular  velocity  vector  of  no-roll  reference  frame  with  respect  to  the 
inertial  frame 
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®D// 

G>D/p 


a> 


B/I 


C/I 


a 


c/i 


nip 

nic 


nip 


a 


D/I 


a 


p/I 


Fp 

Gp/i 


F, 

XaJa,Za 
Y  Y  7 


angular  velocity  vector  of  disk  body  with  respect  to  the  inertial  frame 

angular  velocity  vector  of  disk  body  with  respect  to  the  projectile  body 
frame 

angular  velocity  vector  of  projectile  body  with  respect  to  the  inertial  frame 

translational  velocity  vector  of  two-body  system  with  respect  to  the 
inertial  frame 

translational  acceleration  vector  of  two-body  system  with  respect  to  the 
inertial  frame 

disk  body  mass 

projectile  body  mass 

two-body  system  mass 

nominal  projectile  mass 

translational  acceleration  vector  of  disk  mass  center  with  respect  to  the 
inertial  frame 

translational  acceleration  vector  of  projectile  body  mass  center  with 
respect  to  the  inertial  frame 

weight  vector  of  two-body  system 

weight  vector  of  projectile  body 

weight  vector  of  disk  body 

reaction  force  vector 

angular  acceleration  vector  of  projectile  body  with  respect  to  the  inertial 
frame 

distance  vector  from  projectile  center  of  mass  to  disk  center  of  mass 
aerodynamic  forces  vector 

aerodynamic  force  vector  components  expressed  in  the  body  frame 
aerodynamic  force  vector  components  expressed  in  the  no-roll  frame 
projectile  weight  vector  components  expressed  in  the  body  frame 
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=  projectile  weight  vector  components  expressed  in  the  no-roll  frame 

Z,  F,  Z  =  total  external  force  components  on  the  projectile  body  expressed  in  the 

no-roll  frame 

=  transformation  matrix  from  the  body  frame  to  the  no-roll  frame 

=  reaction  moment  components  due  to  cross  product  of  distance  vector  from 
projectile  center  of  mass  to  disk  center  of  mass  with  reaction  force  vector 
expressed  in  the  body  frame 

Ljj ,  ,  Nj^  =  reaction  moment  components  due  to  cross  product  of  distance  vector  from 

projectile  center  of  mass  to  disk  center  of  mass  with  reaction  force  vector 
expressed  in  the  fixed  frame 

Hpi,  =  angular  momentum  vector  of  the  projectile  body  with  respect  to  the 

inertial  frame 

H jj/j  =  angular  momentum  vector  of  the  disk  body  with  respect  to  the  inertial 

frame 

M  f,  =  reaction  moment  vector 

M  ^  =  total  aerodynamic  moment  vector 

Hp^  ,Hpy,  Hp^  =  angular  momentum  derivative  vector  components  of  the  projectile  body 

expressed  in  the  body  frame 

H  px,H  py,H  pz=  angular  momentum  derivative  vector  components  of  the  projectile  body 

expressed  in  the  no-roll  frame 

Hdj,  ,  Hp,y ,  Hp,z  =  angular  momentum  derivative  vector  components  of  the  projectile  body 

expressed  in  the  body  frame 

Hdx  ,  Hoy  ’  ^Dz  -  angular  momentum  derivative  vector  components  of  the  projectile  body 
expressed  in  the  no-roll  frame 

L^,M =  total  aerodynamic  moment  vector  components  expressed  in  the  body 
frame 

,  N =  total  aerodynamic  moment  vector  components  expressed  in  the  no-roll 
frame 

Ip  =  mass  moment  of  inertia  of  the  projectile  body  about  its  mass  center  with 

respect  to  the  body  frame 
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h 

L,M,N 


C 

'^zo 

C 

'^zo 

c 

^X2 

c 

*^70 

c 


c 

^zo 

c 

'^ZBl 


c 

'^YPA 

c 

'^DD 


c 


LP 


c 


MQ 


V 

a 

P 

P 


mass  moment  of  inertia  of  the  disk  body  about  its  mass  center  with  respect 
to  the  disk  frame 

nominal  mass  moment  of  inertia  of  the  projectile  body  about  its  mass 
center  with  respect  to  the  body  frame 

total  aerodynamic  moment  vector  components  expressed  in  the  no-roll 
frame 

zero  yaw  axial  force  aerodynamic  coefficient  for  the  projectile  body 

zero  yaw  axial  force  aerodynamic  coefficient  for  the  projectile  body 

yaw  drag  axial  force  aerodynamic  coefficient  for  the  projectile  body 

trim  side  force  aerodynamic  coefficient  for  the  projectile  body 

normal  force  aerodynamic  coefficient  along  axis  for  the  projectile 
body 

trim  vertical  force  aerodynamic  coefficient  for  the  projectile  body 

normal  force  aerodynamic  coefficient  along  axis  for  the  projectile 
body 

Magnus  force  aerodynamic  coefficient  for  the  projectile  body 

roll  moment  aerodynamic  coefficient  due  to  fin  cant  for  the  projectile 
body 

roll  damping  moment  aerodynamic  coefficient  due  to  fin  cant  for  the 
projectile  body 

pitch  damping  moment  aerodynamic  coefficient  due  to  fin  cant  for  the 
projectile  body 
velocity  magnitude 

longitudinal  aerodynamic  angle  of  attack 
lateral  aerodynamic  angle  of  attack 
density  of  air 
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ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 
COPIES 


ABERDEEN  PROVING  GROUND 

2  DIR  USARL 

AMSRL  Cl  LP  (BLDG  305) 
AMSRL  Cl  OK  TP  (BLDG  4600) 


1  COMMANDING  GENERAL 
US  ARMY  MATERIEL  CMD 
AMCRDA  TF 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  INST  FOR  ADVNCD  TCHNLGY 
THE  UNIV  OF  TEXAS  AT  AUSTIN 
3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

1  US  MILITARY  ACADEMY 

MATH  SCI  CTR  EXCELLENCE 
MADN  MATH 
THAYER  HALL 
WEST  POINT  NY  10996-1786 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  D 
DR  D  SMITH 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  CS  IS  R 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  OK  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  CS  IS  T 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


2  DEFENSE  TECHNICAL 
INFORMATION  CENTER 
DTIC  OCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FT  BEL  VOIR  VA  22060-6218 
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